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ABSTRACT 

We present results from an 87-ks Suzaku observation of the canonical low-excitation radio galaxy (LERG) 
NGC 6251. We have previously suggested that LERGs violate conventional AGN unification schemes: they 
may lack an obscuring torus and are likely to accrete in a radiatively inefficient manner, with almost all of the 
energy released by the accretion process being channeled into powerful jets. We model the 0.5-20 keV Suzaku 
spectrum with a single power law of photon index T = 1.82^5" o5» together with two collisionally ionized 
plasma models whose parameters are consistent with the known galaxy- and group-scale thermal emission. 
Our observations confirm that there are no signatures of obscured, accretion-related X-ray emission in NGC 
6251, and we show that the luminosity of any such component must be substantially sub-Eddington in nature. 
Subject headings: galaxies: active - galaxies: jets - galaxies: individual (NGC 6251) - X-rays: galaxies 



1. INTRODUCTION 

The origin of relativistic jets in active galactic nuclei (AGN) 
is a key unsolved problem in extragalactic astrophysics. 
While 90% of all AGN (Seyfert galaxies and radio-quiet 
quasars) show little or no jet emission, the remaining 10% 
(the radio-loud AGN and radio-loud quasars) launch power- 
ful twin jets of particles from their cores. Since jets transport 
a significant fraction of the energy liberated during the accre- 
tion process, sometimes out to ^Mpc distances, understand- 
ing how they are produced is critical to a complete picture of 
accretion and feedback in AGN. 

X-ray observations of the nuclei of radio-loud AGN are par- 
ticularly useful for establishing the connection between the 
accretion flow, black hole, and jet. Since the launch of Chan- 
dra and XMM-Newton, several groups have performed spec- 
troscopic studies of the nuclei of 3CRR radio galaxies (e.g., 
Evans et all 120061: iBalmaverde et al.l 120061: lHardcastle et all 
20071 12009b . These efforts have found a dichotomy in the 
nuclear X-ray spectra of radio-loud AGN, related to the opti- 
cal emission line classifications, which may b e related to their 
differ ing roles of hot and cold gas accretion (Hardcastle et al. 
12007 1 ). or possibly the spi n configuration of the black hole it- 
self dGarofalo et alj|2010b . 

High-excitation radio galaxies (HERGs), those with promi- 
nent emission lines in their optical spectra, have X-ray spectra 
that are consistent with unification models. HERGs have stan- 
dard, cool luminous accretion disks and strong (yet typically 
unbroadened) Fe Ka lines. The nuclear continuum is heavily 
obscur ed in the X-ray by cold gas with columns of ~10 
cm~ 2 (Ev ans et al.l 12006) when the source AGN is viewed 
close to edge-on with respect to the observer, and are largely 
unobscured when oriented closer to the line of sight. These 
sources also have an additional component of unabsorbed 
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emission that is associated with the pc- scale jet. H ERGs tend 
to occupy gas-poor environments (e.g., Kraft et al. 2007|) and 
their host galaxi es lie in the 'g reen valley' of the galaxy color- 
mass diagram (Smolcic2009). 

Low-excitation radio galaxies (LERGs), on the other hand, 
the population of radio-loud AGN with little or no nuclear op- 
tical line emission, seem to lack any of the features required by 
standard AGN unification models. We have previously argued 
that their X-ray emission is dominated by radiation from a 
parsec-scale jet, that they have no signatures of standard, cold 
accretion disks, and that t hey may show no evidence at all fo r 
an obscuring torus (lOgle et al.ll2006t lHardcastle et alJl2009h . 
We have suggested that accretion in LERGs takes place in 
a radiatively inefficient manner, with almost all the avail- 
able energy from accretion being channeled into jets. LERG 
host g alaxies lie in t he red sequence on the color-mass dia- 
gram (Smolcic 2009) and are often asso ciated with group - or 
cluster-scale hot-gas environments (e.g.. lTasse et alJ l2008) 

Constraining the properties of the accretion flow in LERGs 
with deep X-ray observations would provide a strong test of 
our hypothesis that the X-ray emission is dominated by a jet 
with no signatures of cold-gas accretion. However, existing 
observations of LERGs have been restricted to relatively poor 
signal-to-noise Chandra and XMM-Newton data, while the 
>10 keV band where the Compton reflection hump could lie 
is completely inaccessible to those observatories. In order to 
address these issues, we obtained an 87-ks Suzaku observation 
of NGC 6251 (z=0.0247), the X-ray brightest low-excitation 
radio galaxy in the 3CRR catalog. We use the XIS and HXD 
to search for signatures that are related to a standard, cool ac- 
cretion disk that is surrounded by a torus: (1) a ~6.4 keV 
Fe Ka line, (2) heavily absorbed continuu m emission, and 
(3) th e >10 keV Compton reflection hump (George & Fabian 
1991). Detection of any of these features would immediately 
invalidate our model. 

Previous attempts to understand the origin of X-ray emis- 
sion in NGC 6251 have yielded conflicting results, making 
our Suzaku observ ation particularly useful. For example, 
iTurner et all d 19971) (based on ASCA data) and iGliozzi et all 
(2004) (based on XMM-Newton data) claimed a detection of 
a high-equivalent width Fe Ka line, suggestive of a dom- 
inant contribution from a luminous accretion d isk. How- 
ever, reanalysis of the XMM-Newton spectra by Eva ns et alJ 
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(2005) instead showed no evidence of an Fe Ka line, which 
favors a jet-dominated scenario. Furthermore, the double - 
peaked SEP of NGC 6 251 measured by, e.g [Hoi (fl999l) 
IGuainazzi et al.1 (120031). IChiaberge et al.1 (120031) . Evans et alJ 
d2005l) . and iMigliori et alj (1201 ll) led those authors to con- 
clude that the high-energy emission is synchrotron self- 
Compton (SSC) emission from a jet. Finally, analysis of the 
spectral vari ability from a long -term monitoring campaign 
with KXTE dGliozzi et al.l 120081) again suggests that the jet, 
rather than the accretion flow, dominates the X-ray emission. 
Nonetheless, only the combination of high effective area and 
simultaneous soft and hard band offered by Suzaku can re- 
solve once and for all the issues surrounding the Fe Ka line, 
as well as place strong limits on the luminosity of the accre- 
tion flow. 

This paper is organized as follows. Section [2] provides a 
summary of our observations and a description of their re- 
duction. In Section [3] we present the results of our spectral 
fitting to the XIS and HXD data. In Section g] we show that 
an Fe Ka line is not detected, while in Section [5] we discuss 
whether or not a buried AGN is present in NGC 625 1 . We end 
with our conclusions in Section[6] 

2. OBSERVATIONS AND DATA REDUCTION 

We observed NGC 625 1 with Suzaku on 2010 December 02 
(OBSID 705039010) for a nominal exposure of 87 ks. Both 
the X-ray Imaging Spectrometer (XIS) and Hard X-ray Detec- 
tor (HXD) were operated in their normal modes. The source 
was positioned at the nominal aimpoint of the XIS instrument. 
The data were processed using v. 2.5.16.29 of the Suzaku pro- 
cessing pipeline, which includes the latest Charge Transfer In- 
efficiency (CTI) correction applied for the XIS. We used the 
standard cleaned events files, which are screened to remove 
periods during which the satellite passed through the South 
Atlantic Anomaly (SAA), had a pointing direction <5° above 
the Earth, or had Earth day-time elevation angles <20°. We 
describe our analysis of the XIS and HXD data below. 

2.1. XIS 

We used data from the two operational front-illuminated 
(FT) CCDs (XIS0 and XIS3), together with the back- 
illuminated XIS1 detector. The three XIS CCDs were op- 
erated with a frame time of 8 s. For our analysis, we used 
data taken in the 3x3 and 5x5 edit modes. We selected only 
events corresponding to grades 0, 2, 3, 4, and 6, and removed 
hot and flickering pixels with the CLEANS IS tool. 

We extracted the spectrum of NGC 6251 from the XIS 
CCDs using a source-centered circle of radius 260" (250 
physical pixels), with background sampled from an adjacent 
region free from any unrelated sources, as well as the 55 Fe 
calibration sources at the corners of each detector. We gen- 
erated response matrix files (RMFs) for each detector using 
v. 2009-02-28 of the XISRMFGEN software, and ancillary re- 
sponse files (ARFs) using v. 2009-01-08 of the XISSIMARF- 
gen software. 

The net exposure times and count rates for the three CCDs 
are shown in Table Q] We co-added the two FI spectra using 
the ADDASCASPEC program, and grouped the resulting spec- 
trum and that of the XIS 1 detector to a minimum of 50 counts 
per bin in order to use \ 2 statistics. We restricted the energy 
range for our spectral fitting to 0.5-10 keV. 
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Fig. 1 .— Suzaku XIS FI (black), XIS BI (red), and HXD PIN (green) best- 
fitting spectra and residuals. The model fit is the sum of a power law with 
slight excess absorption (consistent with the known dusty disk), and two col- 
lisionally ionized plasma components (Model IV). 

We extracted the source spectrum from the HXD/PIN de- 
tector, using the cleaned PIN events files described above. 
The source was not detected with the GSO instrument. We 
used the HXDPINXBPI script, which generates the PIN non X- 
ray background (NXB) spectrum from the latest 'tuned' time- 
dependent instrumental background event file provided by the 
Suzaku Guest Observer Facility. The script also extracts the 
source and background spectra using a common Good Time 
Interval (GTI) criterion, creates a Cosmic X-ray Background 
(CXB) spectrum, and adds the CXB and NXB spectra to- 
gether. We adopted the default binning criterion of 20 counts 
per channel. NGC 6251 is detected at energies between 16 
and 20 keV, and so we restrict our subsequent spectral analy- 
sis to this range. 

3. SPECTRAL ANALYSIS 

In all spectral fits, we linked the parameters of the model 
components across the Suzaku XIS FI, XIS BI, and HXD PIN 
spectra. We tied the normalizations of the FI and BI XIS 
spectra together, but applied a constant cross-normalization 
factor of 1.16 for the PIN spectrum relative to the XIS as 
described in the Suzaku Data Reduction GuideQ. We per- 
formed our spectral fitting using v. 12.6.0 of the XSPEC spec- 
tral fitting package. All results presented here use a cosmol- 
ogy in which Cl m ,o = 0.3, £Ia.o = 0.7, and Ho = 70 km s 
Mpc" 1 . At the redshift of NGC 6251, z=0.02471, the lumi- 
nosity distance is 267.3 Mpc. Errors quoted are 90 per cent 
confidence for one parameter of interest (i.e., Xmin + 2-7), 
unless otherwise stated. All spectral fits include the Galac- 
tic absorption to NGC 625 1 of A^Gai = 1-69 x 10 20 cm" 2 
(iDickey & Lockmanll 19901) . 

We initially attempted to model the spectrum using a single, 
unabsorbed power law, but this resulted in a poor fit (x 2 =1034 
for 832 dof) and noticeable residuals around 1 keV (Model I) 
that suggested t he need for a th e rmal component. Fo l lowin g 
the methods of I Guainazzi et al.l (120031) . iGliozzi et all (120041) . 
and lEvans et al.l (120051) . we subsequently added a collisionally 
ionized plasma (Apec) component with i ts abundance fixed at 
0.35 of solar (similar to lEvans et alj2005l) . which resulted in a 
substantial improvement to the fit (A% 2 =134 for 2 additional 
parameters), but still gave significant residuals below 1 keV 
(Model II). Next, we added a second Apec component to the 
model, which improved the fit to x 2 =87 1 for 828 dof (Model 



2.2. HXD 



' http://heasarc.gsfc.nasa.gov/docs/suzaku/analysis/abc/ 



The Suzaku X-Ray View of NGC 625 1 



3 



I" 



5.5 



6 

Energy (keV) 



6.5 



FIG. 2. — Suzaku FI (black) and BI (red) data/model ratio in the energy 
range 5-7 keV. Also marked is the position where a neutral Fe Kce line would 
lie (dashed line). 

III). Finally, owing to this model's overprediction of flux be- 
low 0.7 keV, we added neutral intrinsic absorption at the red- 
shift of NGC 6251 to the power law component (Model IV). 
This resulted in a substantial improvement to the fit (Ax 2 =33 
for 1 additional parameter). The best-fitting spectrum and 
model are shown in FigureQ] The addition of additional com- 
ponents, such as heavily obscured emission, failed to provide 
a significant improvement to the fitting statistic, and so we 
adopt Model IV as our best fit in the following discussion. 

The measured intrinsic absorption is Nn=(7 .68~^2i$) x 
10 20 cm~ 2 . The detection of modest absorption at the redshift 
of NGC 625 1 is not surprising, given the known dusty disk in 
the host ga laxy. Indeed, the visual e xtinction of Ay = 0.61 ± 
0.12 mag (Ferr arese & Fordlll999 !) corresponds to a neutral 
hydrogen column density of (1.3 ± 0.3) x 10 21 cm~ 2 , which 
is close to the measured value from the Suzaku spectrum. The 
power law photon index of r=1.82+QQg is consistent with the 
range of values fro m other X-ray observatories presented by 
lEvans et ail (120051) . The 2-10 keV unabsorbed luminosity of 

the power law is (2.78+q;24) x 10 42 er 8 s s_1 ' w hich again is 
within the range of historical values discussed by Evans et al. 
(2005). We also detect two collisionally ionized plasmas, with 
temperatures 0.80 ± 0.15 keV and kT=lA2t° i 5 keV. Their 
temperatures are consistent with the extended emission de- 
tected by XMM- Newton on scales of tens and hundreds o f 
kpc, respectively (Sambruna et al. 2004 lEvans et al.l [20051) . 
Furthermore, the fraction of the integrated l uminosity of the 
beta profile determined by lEvans et al.l (12005) within our 260" 
Suzaku extraction region is approximately 60%. This corre- 
sponds to a bolometric luminosity of 4.4 x 10 41 ergs s~\ 
which is consistent with the luminosity we measure with 
Suzaku. 

4. NO EVIDENCE OF IRON Ka EMISSION 

A key goal of our Suzaku observation was to search for 
Fe Ka emission, which would indicate the presence of reflec- 
tion from an accretion disk, and/or circumnuclear torus. Fig- 
ure|2]shows the XIS FI and BI data/model ratios of NGC 625 1 
in the energy range 5-7 keV. It is immediately evident that no 
Fe Ka line is detected with Suzaku. Nonetheless, we added to 
our best-fitting spectral model a Gaussian line with its energy 
fixed at 6.4 keV and its width fixed at 50 eV (below the instru- 
mental resolution of XIS). This resulted in a small improve- 
ment to the fit (Ax 2 =l-38 for 1 additional parameter), with 
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FIG. 3.— Suzaku XIS FI (black), XIS BI (red), and HXD PIN (green) un- 
folded spectra. The overall model fit shown (blue) consists of an absorbed 
power law (pink) and two thermal components (maroon), plus the direct, re- 
flected, and fluorescent components of MYTorus (gray). We have assumed a 
column density of 10 23 cm -2 (upper panel) and 10 24 cm -2 (lower panel). 
The absorbed power law had a best-fit normalization of 0, but for illustrative 
purposes we show its 3cr upper limit. In both cases, we assumed that there 
were no systematic uncertainties in the PIN background. 

a probability of 24% that the improvement is due to statisti- 
cal fluctuations alone. The lower limit to the normalization of 
the Gaussian is consistent with zero, while the upper limit to 
the equivalent width (with respect to the overall continuum) 
is ~100 eV. We therefore conclude that there is no substan- 
tial evidence for neutral Fe Ka emission, although we caution 
that we cannot rule out its presence entirely. 

5 . IS A BURIED AGN PRESENT IN NGC 625 1 ? 

We next used the entire Suzaku bandpass to search for 
the presence of any signatures related to an accretion flow 
and torus. We assumed that , in addition to the dominant 
power-law component (which Evans et al. 2005 interpreted to 
be associated with a jet), there exists a "buried", highly ob- 
scured AGN. We added to o ur best-fitting spectral des cription 
the MYTorus table models dMurphy & Yaqoobll2009h . a fully 
self-consistent treatment of reflection and absorption from a 
circumnuclear torus. Specifically, MYTorus includes compo- 
nents that represent the (absorbed) zero-order continuum, the 
scattered continuum, and resulting neutral fluorescent lines. 
We assumed (1) the primary (accretion-related) power law has 
a photon index of 1 .7 (consistent w ith the mean phot on index 
of z < 0.1 3CRR sources found by lE vans et al.l2006b : and (2) 
it is attenuated by equatorial column d ensities of 10 23 cm " 2 
(consistent with HERG s measured by Evans et al. 20061 and 
Balm averde et alj|2006l) or 10 24 cm -2 (to test the Compton- 
thick case). The choice of inclination angle is problematic: 
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iJones & Wehrld d2002l) used the jet-counterjet brightness ra- 
tio to conclude that the angle to the line of sight, 6, of the jet 
is ^45°, b ut this is different to th at subtended by the dusty 
disk (76°) dFerrarese & Fordll 19991) . For illustrative purposes, 
we adopt an angle of 90° (i.e., an edge-on geometry) and 
note that our results are essentially unaffected by the choice 
of angle until 9 < 60 ° , after which the reprocessor adopted 
by Murphy & Yaqoob (2009) no longer intersects our line of 
sight. 

The XIS and HXD spectra and models for the 10 23 cm~ 2 
and 10 24 cm -2 cases are shown in Figure [3] The best- 
fitting accretion-related power law normalizations were con- 
sistent with zero, with 3 a upper limits to the 2-10 keV un- 
absorbed luminosity of 7 x 10 41 ergs s _1 and 4 x 10 42 
ergs s _1 , res pectively. Assuming a black hole mass of 
6 x 10 8 M dFerrarese & Fordl [T999). we find that any ob- 
scured, accretion-related emission is highly sub-Eddington 
(L x /L Edd < 10~ 5 for /Vh=10 23 cm" 2 ; L x /L Edd < 5 x 
10" 5 for 7Y H =10 24 cm- 2 ). The upper limit to the 2-10 keV 
luminosity in both cases is consistent with the relationship be- 
tween X-ray luminos ity and 15/im luminosity established by 
lHardcastle et"ail d2009t . 

As an additional test, we considered an extreme Compton 
thick obscurer {Nn=10 25 cm -2 ). Again, this failed to pro- 
vide a significant improvement to the fit. The upper limit to 
the 2-10 keV unabsorbed luminosity in this case is 2 x 10 44 
ergs s _1 . However, we note that (1) the peak of the Compton 
reflection hump lies outside of the usable energy range of the 
PIN data, meaning that the principal feature used to determine 
the strength of the accretion-related emission is inaccessible, 
and (2) a lu minosity of 2 x 10 44 erg s s _1 would be a 6a out- 
lier from the Hardcastle et al. (2009) L x -Li^^ m relationship. 
Observations with NuSTAR and Astro-H would provide fur- 
ther constraints. 

Finally, given the poor signal-to-noise of our PIN data, it 
is important to address the known systematic uncertainties of 
the Su zaku background, whi ch are of order 3% at the 1-sigma 
level (Fukazawaet aD 120091) . In the above analysis, we no- 
ticed several points in the PIN band that appear to lie above 
the model. However, increasing the background by 3% (set 
using the cornorm parameter in Xspec) removes any discrep- 



ancy, ruling out any evidence of a hard excess given the PIN 
systematics. 

6. CONCLUSIONS 

We have presented results from a new 87-ks Suzaku 
observation of the prototypical low-excitation radio galaxy 
NGC 625 1 . We have shown the following: 

1 . We model the 0.5-20 keV Suzaku spectrum with a sin- 
gle power law of photon index T = 1.82+qq 4 , together 
with two collisionally ionized plasma models whose 
parameters are consistent with the known galaxy- and 
group-scale thermal emission. 

2. There is no significant evidence in the X-ray for ob- 
scured, accretion-related emission. However, we can- 
not rule out accretion-flow luminosities of <7 x 10 41 
ergs s~ 4 ergs s _1 and 4 x 10 42 ergs s _1 if they are 
obscured by columns of 10 23 cm~ 2 or 10 24 cm , re- 
spectively. These luminosities are consistent with those 
predicted by the L x -L^ ,, m relationship established by 
Hardcastle et al] (120091) . Both luminosities are highly 
sub-Eddington (ix/^Edd ~ 10~ 5 ). 

3. NGC 6251 is thus consistent with our predictions for 
low-excitation radio galaxies: they have radiatively in- 
efficient accretion flows, show no evidence in the X- 
ray for circumnuclear tori, and their X-ray emission is 
likely to be jet-dominated in nature. 
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TABLE 1 
Suzaku Observation Log 



Date of observation 


Obs. ID 


Instrument 


Screened exposure time 


Net nuclear count rate (s 1 ) 


2010 Dec 02 


705039010 


XIS0 


87 ks 


0.15±0.02 






XIS1 


87 ks 


0.21±0.02 






XIS3 


87 ks 


0.16±0.02 






HXD/PIN 


74 ks 


(6.46 ±1.18) x 10~ 2 



TABLE 2 

Best-fitting spectral parameters 



Model 


Description 


Nh (cm- 2 ) 


Power Law 


Thermal 1 


Thermal 2 




X 2 /dof 


I 


PL 




r= 1.87 ±0.02 
norm=(6.79±0.09)x 10~ 4 








1034/832 


II 


PL+TH1 




r=1.75±0.03 
norm=(5.62±0.22)x 10~ 4 


fcT=1.2lto.i keV 
Z=0.35 (f) 

norm=(3.17±0.65)xl0 -4 






900/830 


III 


PL+TH1+TH2 




i -i.oz_ 09 

norm=(3.76+° x 10~ 4 


/cT=0.89±0.10keV 
Z=0.35 (f) 


feT=2.63+ ™ keV 
Z=0.35 (f) 




871/828 








norm=(2.04t^; 4 ^) X 10" 4 


norm=(7.86t2 2 g) x 


10" 4 




IV 


ABS(PL)+TH1+TH2 


(7.68+^) X 10™ 


1 -1.82_ 05 

norm=(6.17t ; 5 5) x 10~ 4 


fcT=0.80 ± 0.15 keV 
Z=0.35 (f) 

norm=(1.27j: ;g|) X 10 -4 


feT=1.42+ ; 4 o keV 
Z=0.35 (f) 

norm=(2.22t » 4 ) X 


10" 4 


838/827 



a Col. (1): Model number. Col. (2): Description of spectrum (Abs=Neutral absorption, PL=Power Law, TH=Collisionally ionized plasma model). Col. (3): 
Intrinsic neutral hydrogen column density. Galactic absorption has also been applied. Col. (4): Power Law parameters. Normalization is quoted at 1 keV in units 
of ph keV -1 cm~ 2 s _1 . Col. (5): Parameters of first collisionally ionized plasma model. Col. (6): Parameters of second collisionally ionized plasma model. 
Col. (7): Value of x 2 and de grees of freedom. 



